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We study the dynamics of vortices in ideal and weakly interacting Bose-Einstein condensates using 
a Ritz minimization method to solve the two-dimensional Gross-Pitaevskii equation. For different 
initial vortex configurations we calculate the trajectories of the vortices. We find conditions under 
which a vortex-antivortex pair annihilates and is created again. For the case of three vortices we show 
that at certain times two additional vortices may be created, which move through the condensate 
and annihilate each other again. For a noninteracting condensate this process is periodic, whereas for 
small interactions the essential features persist, but the periodicity is lost. The results are compared 
to exact numerical solutions of the Gross-Pitaevskii equation confirming our analytical findings. 



PACS numbers: 03.75.Kk, 03.75.Lm, 05.30.Jp, 32.80.Pj 



I. INTRODUCTION 

Quantized vortices play an important role in verify- 
ing the superfluid properties of quantum liquids such as 
Bose-Einstein condensates (BECs) or degenerate Fermi 
gases. In weakly interacting alkali gases condensate 
states containing a single vortex line were first created us- 
ing Raman transition phase- imprinting methods By 
rotating the system with a laser spoon 0, 0] , vortex lat- 
tices containing more than 100 vortices have been created 
UGL and by using topological phase engineering meth- 
ods @ it is even possible to create multiply charged vor- 
tices. It is expected that more complicated vortex clus- 
ters can be created in the future, e.g., with the further 
development of phase- imprinting methods. Such states 
would enable various opportunities, ranging from inves- 
tigating the properties of randompolynomials Q to using 
vortices in quantum memories [8|]. All of these develop- 
ments stir a great interest in the study of states with 
several vortices. 

Recently, there were a number of investigations on the 
properties of vortices in BECs. For three-dimensional 
condensates, several studies on the dynamics of vortex 
lines have been done [9|, [1CJ . The generation and dynam- 
ics of vortices in a toroidal condensate have been inves- 
tigated in Ref. [Tlj] , whereas detailed numerical studies 
of the optical generation of vortices in pancake-shaped 
condensates have been carried out in [12j . The manip- 
ulation of vortices such as charge conversion by exter- 
nal potentials has been discussed in Ref. [H|. Further 
numerical studies revealed that for condensates with a 
strong nonlinearity there exist several configurations of 
vortices which are stable [lH [TEI . flfjj . Analytical ex- 
pressions for the angular momentum and the energy of a 
vortex-antivortex configuration in a BEC have been ob- 
tained [l~7j using the Thomas-Fermi approximation. For 
strongly nonlinear condensates analytical solutions were 



derived by splitting the wave function into a region close 
to the vortex core and one far away from the vortex, 
where the hydrodynamic properties of the condensate are 
an essential feature [lj| H3| • In contrast, here we con- 
centrate on the dynamical properties of vortex configura- 
tions in weakly interacting condensates, and solutions for 
the whole spatial regime are obtained. In a noninteract- 
ing condensate the vortices behave similar to those cre- 
ated in an optical beam using holograms [U [22j , where 
the time evolution in the BEC corresponds to the spatial 
evolution of the laser beam. As we will show, the inter- 
action between the BEC atoms changes the behavior of 
the vortex dynamics considerably. 

In this paper, we make use of the Gross-Pitaevskii 
equation (CPE), also known as the nonlinear Schrodinger 
equation, which is known to be a valid description of the 
mean-field dynamics of a BEC at zero temperature. We 
consider a harmonic trap with tight confinement along 
one direction, such that the condensate is effectively two 
dimensional. For the case of an ideal, i.e., noninteracting 
BEC, the dynamics of the vortices is solved analytically 
yielding the essential features of the time evolution. For 
a small interaction within the condensate, we make use 
of the Ritz method in order to get analytical estimates 
of the dynamics. These estimates are compared to exact 
numerical solutions of the CPE using the time-splitting 
spectral method (TSSP), which is explicit, uncondition- 
ally stable and spectrally accurate in space. Details of 
the numerical method are described in Refs. [23l. |24|. l2a|. 

This paper is organized as follows. In Sec.|TT]we intro- 
duce the model under investigation and define the general 
initial states of the vortex configurations. In Sec. IIIII we 
discuss the dynamics of vortices in an ideal condensate 
as a background for the results of Sec. IIV1 where a de- 
tailed investigation of the dynamics of a single vortex, a 
vortex pair, a vortex dipole, and a vortex tripole are pre- 
sented. We conclude in Sec. [Vj An Appendix contains 
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some more details on the evolution of a vortex tripole for 
the noninteracting BEC. 



II. THE MODEL 

In this work, we consider a Bose-Einstein conden- 
sate (BEC) in a radially symmetric trap Vt(x, y, z) = 
^mb[u>(x 2 + y 2 ) + uo z z 2 ] with u> z ^> to the trap frequen- 
cies in axial and radial direction, respectively, and raj, 
the mass of the BEC atoms. We assume a tight con- 
finement in axial direction such that Tilo z 3> k^T, where 
&b is Boltzmann's constant and T is the temperature of 
the BEC, as well as Tiw z 3> gn^, with no the density 
of the BEC in the center of the trap and g the interac- 
tion strength within the BEC, given by g = 4:Trh 2 a s /mb : 
with a s the s-wave scattering length. For tempera- 
tures well below the critical temperature of the BEC 
and the Berezinskii-Kostcrlitz-Thouless transition tem- 
perature [2rl I27I , a regime which is in reach of current 
experiments [28|, l29| , phase fluctuations occur on scales 
which are typically larger than the size of the condensate 
and the dynamics of the BEC is well described by the 
dimensionless 2D Gross-Pitaevskii equation 23] 
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Here, tp = ip(x,y,t) is the normalized wave function 
of the condensate with J \tp(x, y)\ 2 dxdy — 1, and (3 = 
2Na s yj2Tiuj z /uo I clq characterizes the interatomic inter- 
action, defined in terms of the total number of particles 
N in the condensate. The above dimensionless quanti- 
ties are obtained by scaling the length by the harmonic 
oscillator length ao = y/fi/m^uj, the time by and 
the energy by Tko. 

We study the dynamics of n vortices with topological 
charge qj = ±1 (j = 1, 2, n), which are initially placed 
at positions Tj — (xj,yj). For this purpose, we first need 
to calculate the ground state tp gs (r) of the GPE and the 
state ijjq (r) with a single vortex in the center of the trap. 
From this, we extract the function p q (r) — i/' 9 (r)/V'gs(r) I 
which describes a vortex in the BEC background. The 
initial state for the vortices is then approximately given 
by 



i>(r,t = 0) =aVg 8 (r) JJpg 3 (r-rj) 



(2) 



where a is chosen such that the initial state is normal- 
ized to 1. Unless otherwise stated, this normalization 
constant will be dropped in the following. The above ap- 
proximation holds for vortices which are not too close to 
the edge of the condensate. A physical realization of such 
states can be achieved either by stirring the condensate 
@; HI or, in a more controlled way, by phase imprint- 
ing methods @ or by using light with orbital angular 
momentum 



III. VORTICES IN A NONINTERACTING BEC 

To get an insight into the dynamics of vortices and as 
a background against the results found for the weakly 
interacting condensate we first focus on the noninter- 
acting case j3 = 0, where the GPE simplifies to a 
two-dimensional harmonic oscillator. The initial state 
tp(r, t = 0) is expanded in terms of the solutions 



tp n , m (r,t) = 



-iE„, m t-(a; 2 +y 2 )/2 



\/2 n + m nlmlir 



H n (x)H m (y) (3) 



of the harmonic oscillator, with the energies E n>m = 1 + 
n + m, H n are Hermite polynomials, and integer numbers 
n,m> 0. This also gives the time evolution of the vortex 
state. The trajectories of the vortices are calculated by 
finding the zeros of the wave function and checking if 
at these points the condensate has a nonzero winding 
number. 

For a single vortex of topological charge q = +1 ini- 
tially located at (x\ , 0) the expansion gives the wave func- 
tion 
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2it\ [x + iy - e u xi] 



<i>v{x,y,t) = exp 

The position of the vortex evolves in time on an exact 
circular trajectory described by x v (t) — xi cos(i) and 
y v (t) = xi sin(<). 

For a vortex pair, that is two vortices with identical 
topological charge q = +1 initially located at (xj,yj), 
where j — 1,2, the expansion into the solutions of the 
harmonic oscillator shows that the two vortices move in- 
dependently from each other on trajectories x vp j(t) = 
Xj cos(i) — yj sin(i) and y V p,j(t) = Xj sin(i) + yj cos(i), 
which are exactly the same as for a single vortex. Due 
to the conservation of the total topological charge Q = 
J2j Qj there are always at least two vortices present in 
the condensate, since a doubly charged vortex is unstable 
[30| . From the absence of additional vortices we conclude 
that the total energy of the system is not large enough 
to allow the spontaneous creation of a vortex-antivortex 
pair. 

The dynamics changes considerably for a vortex dipole. 
Assuming the q = +1 vortex is initially located at 
(xq,0) and the q = — 1 vortex at (— xo,0), the tra- 
jectories are given by y v d(t) — sm(t)(xQ — 1)/xq and 
x v d{t) = ±\/xo — 2/vd W- Examples of these trajectories 
are shown in Fig. OJa). For xo > l/v2 these expressions 
are always real numbers indicating that the two vortices 
do not annihilate each other. For xq = 1 the vortices 
remain even stationary at their initial positions. The sit- 
uation changes for xq < \j\[2. In this case, the vortices 
collide with each other at a time t a given by the first 
possible solution of 



|sin(t)| = 



Xq 



(5) 
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FIG. 1: (a) Trajectories for a vortex dipole for different 
initial positions (ixo, 0). From inside to outside xq — 
0.4, 0.6, 0.8, 1, 1.2. For x = 0.4, 0.6, the two vortices annihi- 
late each other at certain times and reappear again, indicated 
by the closed lines. For all other cases shown, there are always 
two vortices present apart from times where t = (2n + l)7r/2, 
with n an integer number, (b) Trajectories for the more gen- 
eral initial condition ri = (0.9,0.1), r-2 = (—1.1,0). 



They annihilate each other and only reappear again at 
time t r , where the second possible solution of Eq. ([5]) in 
the half-period t <E [0, n] is taken. Between those two 
times, there are no vortices present in the condensate. 
We also did calculations with more general initial con- 
ditions for the two vortices, located at rj = (xj,yj), as 
shown in Fig.[TJb). The whole dynamics is still periodic 
with a period of 27r, and after a time of t = it the state 
is given by (f> v d{x , y , ir) — — </> v d(— x, — y, 0), which means 
the initial state is, up to an unimportant global phase, in- 
flected at the origin. For the symmetric stable case with 
Xq = 1 this might be surprising at the first glance, but 
is explained by the fact that the (dimensionless) current 
density j(x, y, t) — — i (i/}*\/tp — ipVip*) /2 stops at times 
t = (2n + l)w/2, n = 1, 2, allowing the two stationary 
vortices to flip their signs. 

The dynamics of a vortex tripole with two vortices 
of topological charge q = +1 at locations (xo,0) and 
(— xq, 0) and one of charge q = — 1 at (0, 0) is given by 



<Pvt(x,y,t) 



[x 3 + iyx 2 + (y 2 - e 2it (x 2 - 2) - 2) x (6) 
+ i y ( ? y 2 + e 2it (x 2 + 2)-2)]. 



For the special case xo = \/2, we find as the ze- 
ros of the wave function (0,0), (-y/2,0), (y/2,0), 
(a/2 - 4cos(2i) sin(2t),-- v /2 - 4cos(2i) cos(2f)), and 
(-a/2 - 4cos(2i) sin(2i), ^2 - 4cos(2i) cos(2t)). Al- 
though the initial state only contains three vortices, 
during the evolution additional vortices are created and 
annihilated again, so that at certain times there is a 
maximum of five vortices present in the condensate. 
The whole evolution is periodic with a period of n, a 
property which will be lost for small interactions. For 
more details of the vortex dynamics in the ideal case we 
refer to the Appendix. 



IV. VORTICES IN A WEAKLY INTERACTING 
BEC 

The situation gets more complicated when a finite in- 
teraction (3 of the condensate is taken into account, and 
exact analytical solutions are not known. In order to cal- 
culate the dynamics of the vortices we therefore proceed 
using the Ritz minimization method 31\. For small in- 
teractions (3 < 1 we assume that the Gaussian shape of 
the condensate is not changed, but only broadened. For 
the solutions of the Gross-Pitaevskii equation we make 
the ansatz 



1pn,m,p(*,t) 



\/2"+ m n!m! 



7T<X 



H n (x/a)H m (y/a), (7) 



where the functions are normalized to 1 and tr is a con- 
stant which takes the broadening into account. This con- 
stant is derived by minimizing the Gross-Pitaevskii en- 
ergy functional [3fJ 



EM = N 



dr||V^| 2 



\r 2 M 2 



(8) 



with respect to tr, where tp = "0o,o,/3 ( r ; 0) and the en- 
ergy is given in units of Huj. The minimum is found for 
a 2 = \J ((3 + 2ir)/2it. The constants fj, n ,m are derived by 
putting the ansatz into the GPE (jXJ) , multiplying by ip* 
and integrating over space, i.e., 



dr-|Vi/> 



■0\tf>. 



'n,m,0 | 



(9) 

The wave functions ij) n ,m,p ( r ! are then used to expand 
the vortex state, which is given by the functions describ- 
ing the vortices multiplied by ipo,o,/3 and a normaliza- 
tion factor. This factor is, as in the preceding sections, 
dropped unless otherwise stated. 

We have further assumed that the shape of the vor- 
tices does not change due to the increased interaction. 
Normally, in a strongly nonlinear condensate the size of 
a vortex is given by the coherence length £, for which 
we find (in scaled units) £ 2 = y^7r/4/3. However, for our 
trapped weakly interacting condensate this length scale 
is no longer useful, which gets especially apparent for a 
vanishing interaction j3 — > 0, leading to £ — * oo. Instead 
by using numerical calculations we find that the vortices 
are still well described by the functions p q derived for 
(3 = 0. This is illustrated in Figs. [Ha) and [21c), where 
the trajectories for a vortex dipole with the analytical 
p q for (3 = are compared to those with numerically 
found p q for (3 = 1. The difference between the trajec- 
tories is negligible. In Figs. [2Kb) and [2{d) we compare 
the vortex trajectories for the factor p q found for (3=1 
and (3 = 2. Again, the trajectories are almost identical. 
This shows that a small perturbation in the shape of the 
vortices leaves the trajectories essentially unaffected as 
long as an initial condition described by Eq. ^ is used. 
We note that Eq. @ only describes initially factorized 
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FIG. 2: (Color online) Comparison of vortex trajectories for 
different initial states for a vortex dipole with initial xo = 1 
((a) and (b)) and x = 0.5 ((c) and (d)). For x < the 
trajectories for a vortex dipole with initial p q for (3 = 1 is 
shown, whereas for x > the trajectories with initial p q for 
/3 = ((a) and (c)) and (3 = 2 ((b) and (d)) are plotted. All 
trajectories are mirror symmetric to the axis where x = 
(dashed line). The vortices were evolved for a time interval 
[0, 20] and a condensate with f3 — 1. 



vortices, which are not necessarily preferred as the ini- 
tial states with a given number of vortices. Alternative 
states can be considered, and it has been shown that 
the dynamics for globally linked, that means nonfactor- 
ized, vortices can behave significantly different from the 
dynamics of factorized ones [32]. However, a full char- 
acterization of such states is beyond the scope of this 
paper and we therefore restrict our considerations to ini- 
tial states described by Eq. 




FIG. 3: (Color online) (a) Precession frequency of a single 
vortex versus the interaction strength (3. Stars show ui p for a 
different initial distance of the vortex from the center, namely 
(from bottom to top) xo = 0.1,0.5,1.0,1.5,2.0. The circles 
show the analytical result given in Eq. (|13p . lines are guides 
to the eye. (b) The coefficients c(xo) for different xo, see text. 



the precession frequency is given by 

(/3 + 8tt) 



(13) 



A comparison with numerical results shown in Fig. [3Ja) 
illustrates that the analytical result describes the trend 
of a decreasing precession frequency well, however, the 
numerical results show a clear dependence on the dis- 
tance xo of the vortex to the center of the trap, which 
is missing in the analytical formula. For small inter- 
action p the behavior of the precession frequency can 
be assumed to be linear with /3 and follow the curve 
Lu p = 1 + c(x Q )/3. For the analytical formula we get 
Cana(^o) = — l/87rw —0.04, whereas the results from the 
numerics are shown in Fig. [3jb). A trend towards the 
analytical value is visible with increasing xq. However, 
measuring the frequency for xq > 2 will become increas- 
ingly difficult due to the dilute condensate density for 
large distances. 



A. Single vortex 



B. Vortex pair 



For a single vortex initially located at (xq,0) the ex- 
pansion into the solutions ipn,m,f3 yields 



<p v (x,y,t) 



i(7g + 16ir)t 



2 , 2 
c +y 

5^3 



x — e iV^-/wF^ X q + iy], 



which immediately leads to the trajectory 

( (13 + 8TT)t 

x v (t) = Xo cos ' 
y v (t) = x sm 



4V27TV/3 + 2tt 
08 + 8Tr)t 



(10) 

(11) 
(12) 



The vortex again moves on an exact circular line, however 
the time it needs to complete one circle is increased com- 
pared to the interaction-free case. To be more specific, 



Let us assume that at time t = the two vortices with 
topological charge qj = +1 are located at (— xo,0) and 
(xo, 0). The time evolution of the BEC wave function is 
given by 



i(137ff + 384ir)t 
64\/27r\/7?+2~7r 



i(41|3 + 25l,- if 



2ie 64 ^ 7 v 



<yx 



e 64v^vTJT2lF X Q - y 



(14) 



We see that the factor describing the two vortices in gen- 
eral cannot be factorized as in the interaction-free case, 
indicating that the two vortices influence each other. 
This also becomes evident when investigating the trajec- 
tories of the two vortices. They no longer move on exact 
circular lines, but rather on deformed ones as shown in 
Fig. 01 The deviation from the exact circular line gets 
larger for higher interaction strength (3. 
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FIG. 4: Trajectories of a vortex pair deduced from the analyt- 
ical formula Eq. (|14|l for /3 = 2, xo = 1.2, and a time interval 
of t 6 [0,20]. The trajectories are no longer exact circular 
lines, indicating that both vortices influence each other. 




FIG. 5: (Color online) Trajectories of the two vortices in a 
vortex dipole for (a) xq — 0.6, (b) Xo — 0.8, (c) xo = 1, 
and (d) xo = 1.2. Analytical results calculated with Eq. (|15p 
are only shown for x < 0, whereas numerical results are only 
shown for x > 0. Both results are mirror symmetric with 
respect to the axis where x — (dashed line). The interac- 
tion is chosen as f3 = 1, and the time interval for which the 
trajectories are shown is given by < t < 20. 



C. Vortex dipole 

For symmetric initial conditions, i.e., a negative vortex 
initially at position (— xq, 0) and a positive one at position 
(xq, 0), the wave function for a vortex dipole is expanded 
in terms of ipn,m,(3 yielding 



0vd(r,t) 



3i(115g + 256x)t 



i(233g + 512ir)t 
2i e 64V^V75T21F XQ y 



i(249/3 + 6407r)t 

_|_ e 6 4v^rvTJ+2iF (fj _ a; )(cr 

i(13ff + 24ir)t 




(15) 



As for the vortex pair the dynamics of the BEC is 
no longer periodic due to the interaction, but acquires 
a more complicated time dependence. As an example, 
the trajectories for f3 = 1, a time interval of t € [0,20], 
and several initial positions (±xo,0) are shown in Fig. [5] 
The simple trajectories from the noninteracting case are 
changed to complicated structures, which have lost their 
periodicity. 

Comparison between analytical results using Eq. (|15|) 
and numerics shows that especially for small times and 
intermediate distances 2xq rs 2 both trajectories agree 
quite well, however for certain times the analytical re- 
sults predict vortices at positions where there should be 
no vortices according to numerics. The differences be- 
tween analytics and numerics get larger for increasing 
interaction (3. Our numerical calculations furthermore 
suggest that even for distances 2x$ > 2 there exist times 
for which the two vortices annihilate each other and reap- 
pear again, so, for example, at 2x = 3 and a time around 
t = 18. This is in contrast to the noninteracting case, 
where the evolution was strictly periodic and an annihi- 
lation of the vortices not possible if their initial distance 
was 2x > 2/a/2. 



D. Vortex tripole 

For a vortex tripole as introduced in Sec. IHII thc wave 
function describing the time evolution is given by 



20i^/27t g £ +x 2 + 2 
871igt SPT^ . - 

(f> vt (x,y,t) —e i28v^/7?+2¥ 6 256^^3+^ ( x + jy) _ 2ix{x - iy)y) 



N (16) 

3i(361/3+10247r)t i(647/il + 2048j)t \ 

_l_ e 256v^v^+2¥ (30-2(2; + iy) - 2 (x 1, + iy 3 )) - 2e i28^2ifv7J+2¥ (2a 2 (x + iy) - (x - iy)xl) . 



As for the noninteracting case, there always exists a vor- able to identify any initial condition where more than the 
tex in the center of the condensate, however we were not 
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FIG. 6: (Color online) (a) Number of vortices N present in 
the BEC vs time, with initially xo = 1.5. (b) Average number 
of vortices (N) present in the BEC during the time interval 
[0,20] for different initial xq. The solid line shows the results 
derived from the analytical expression, Eq. (|16[1 . whereas the 
plus signs show numerical results. In all cases we have chosen 
(3 = 1, apart from the dashed line, which shows the results 
for (3 = 0. 



central vortex are stationary. This is consistent with the 
results reported in Ref. [16]. There it was shown that a 
stable vortex tripole, i.e., a configuration of exactly three 
stationary vortices, only exists for interactions /? > 1, 
where our ansatz is no longer valid. 

The zeros of the wave function Eq. (|16p are found nu- 
merically for different initial positions of the vortices. For 
short times t < 10, the time evolution of the trajectories 
is quasiperiodic, however this quasiperiodicity is more 
and more washed out for longer times. This gets also ap- 
parent when the number of vortices during time is con- 
sidered. As shown in Fig. Oa), the number of vortices 
oscillates for t < 10 between 3 and 5 for an initial con- 
figuration with xq = 1.5. For later times, however, there 
are intervals in which the condensate only exhibits one 
vortex. 

In Fig. [|3b) we show the average number of vortices 
(TV) in the condensate during the time interval [0, 20] ver- 
sus the initial positions xo. The results derived from the 
numerical evolution of the GPE agree quite well with the 
ones from the analytical formula Eq. (fTfj|) . We observe 
that for xo < 1.5 the average number of vortices (N) 
increases with increasing xq, and reaches a pronounced 
maximum at xq ~ 1.5. For larger xq the average vortex 
number decreases again to reach a constant value of 3 for 
xo > 2. This indicates that for such large distances the 
energy within the condensate is too low to spontaneously 
create additional vortices. This behavior is qualitatively 
similar to the noninteracting case as also indicated in 
Fig. [6]^b) , where the interacting case tends to exhibit a 
higher average vortex number. 



V. CONCLUSION 

In the present paper we have investigated the dynam- 
ics of vortices in two-dimensional Bose-Einstein conden- 



sates. We have solved the Gross-Pitaevskii equation an- 
alytically for ideal condensates and have used the Ritz 
minimization method in order to calculate the dynamics 
for small interactions. The latter results were compared 
to exact numerical solutions of the condensate dynamics. 

For an ideal condensate we have shown that two vor- 
tices with the same topological charge \q\ = 1 do not 
influence each other and behave like two independent, 
single vortices. This changes as soon as the interaction 
within the condensate is taken into account. The two vor- 
tices no longer move on exact circular lines, but rather 
on a distorted circular path. 

For the case of a vortex dipole, that is two vortices with 
opposite topological charge in the condensate, we found 
that for certain initial conditions in an ideal condensate 
the two vortices will collide, thereby annihilating each 
other, and reappear again, whereas they are always sep- 
arated if the initial distance is large enough. The trajec- 
tories change considerably in an interacting condensate, 
and a large initial distance will no longer guarantee that 
the two vortices do not annihilate each other at some 
times. 

We also investigated the case where initially there are 
three vortices in the BEC, a so-called vortex tripole. In 
contrast to the case of only two vortices present at the 
beginning we found that for an ideal condensate during 
the time evolution additional vortices were created and 
annihilated, allowing for a maximum of five vortices in 
the condensate. Our numerical results showed that this 
behavior also persisted in an interacting condensate. 
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APPENDIX A: DETAILS OF THE VORTEX 
TRIPOLE DYNAMICS FOR THE IDEAL BEC 

The expressions given for the vortex trajectories in 
Sec. IIIII do not include information on the topological 
charge of the respective vortices. Closer investigations 
show that the vortices described by a single trajectory 
can indeed flip their signs. To be more precise, we as- 
sume that at t = the central vortex has a negative 
charge. At times t < 7r/6, there are only three vortices 
present in the BEC, which do not move. For t = 7r/6, 
two new vortices of equal charge q = — 1 arise in the cen- 
ter of the BEC flipping the charge of the central vortex 
in order to keep the total topological charge Q constant. 
The two new vortices start to move out of the center to- 
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FIG. 7: (Color online) (a) Trajectories of the two additional 
vortices (black lines) for the initial state of a vortex tripole 
with xo — \/2, details see text. The three stationary vortices 
are indicated by plus signs, (b) Density plot for t = 0.2tt. 
Plus signs indicate the position of the five vortices. 



wards the two stationary vortices. When the two new 
vortices cross the two stationary ones at t — 7r/4, their 
charges are flipped as well, such that the moving vor- 
tices have now a positive charge and the stationary ones 
a negative one. After the flip, the nonstationary vortices 
move around the three stationary vortices as indicated in 
Fig. EJa), cross the two vortices outside the center again 



and flipping their charge at t = 3tt/4, and finally anni- 
hilate at the center at a time t = 5ir/6, such that the 
original configuration for t — is achieved again. This 
is repeated periodically, with a period of 7r. Investigat- 
ing the phase of the condensate at the times t = ir/4 and 
37r/4, where the moving vortices cross the stationary ones 
at (±V2, 0), shows that for these times there is only one 
vortex present in the condensate, namely the central one. 
A density plot of the system for a time t = 0.2ir is shown 
in Fig.[7Jb). The additional vortices are created in a re- 
gion of low density, where the energetic cost is lowest. 
However, this low density makes the measurement of the 
vortices more difficult. 

For the case xq ^ V2 we find a similar behavior of the 
vortices. For xo > 2 there will always be three vortices 
in the condensate. For smaller values, additional vortices 
can be created, and for xq < \/2 we find that for certain 
time intervals (i.e., not only for single points in time) 
there is only one vortex present in the condensate. Due 
to the conservation of the total topological charge this 
is the minimum number of vortices present. Calculating 
the roots of Eq. ^ shows that the maximum number of 
vortices is five. 
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